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Abstract. We consider a Hermite-type basis for which the eigen- 
value problem associated to the operator Ha,b '■= B{—dl) + Ax^ 
acting on L^(M; C^) becomes a three-terms recurrence. Here A and 
B are 2x2 constant positive definite matrices. Our main result 
provides an explicit characterization of the eigenvectors of Ha,b 
that lie in the span of the first four elements of this basis when 
AB ^ BA. 



1. Introduction 

It is well known that the spectrum of the harmonic oscillator Hamil- 
tonian 

Ha ■■= -dl + a^x^, a > 0, 
acting on L^(R) consists of the non-degenerate eigenvalues 
{a{2n+ l)}.^o with corresponding normalized eigenfunctions 

(1) = , r- ' 

where hn{x) := (— l)"e^^(9"[e~^^] is the n-th Hermite polynomial. The 
present paper is devoted to studying the spectrum of a matrix version 
of Ha, the operator 

Ha,b ■.= B{-dl) + Ax\ 

acting on L^(M;C^), where A and B are two 2x2 constant positive 
definite matrices. 

In contrast to the scalar situation, the spectral analysis of Ha,b is 
far more involved due to the non-commutativity of the coefficients. If 
AB = BA, it is not difficult to find the eigenvalues and eigenfunctions 
of Hab from those of Ha- On the other hand, when AB ^ BA, the 
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eigenvalues and eigenf unctions of Ha,b are connected to those of Ha in 
a highly non-trivial manner (see Theorem |21 below). 

Our recent interest in describing spectral properties of operators such 
as Hyi B arises from two sources. In a series of recent works, Parmeg- 
giani and Wakayama, cf. jl], and jH], characterize the spectrum of 
the operator K^^^ := A{-dl + x^) + B{2xd^ + x'^) acting on L'^{R; C^), 

assuming that A is definite positive and B = —5*. Although the two 
operators are related, it does not seem possible to obtain the eigenval- 
ues of Ha^b from those of ^. In j2] and the eigenf unctions of ^ 
are found in terms of a twisted Hermite-type basis of L^(M; C^). In this 
basis the eigenvalue problem associated to ^ becomes a three-term 
recurrence. The strategy presented below for analyzing the spectrum 
of Ha.b will be similar. 

Our second motivation is heuristic. It is known that the scalar har- 
monic oscillator. Ha, achieves the optimal value for the constant in 
the Lieb-Thirring inequalities with power o" > 3/2, cf. [2]. It would 
be of great interest finding Hamiltonians with similar properties for 
Lieb-Thirring-type inequalities for magnetic Schrodinger and Pauli op- 
erators, cf. PP, [21 and [3J. Due to their close connection with the 
harmonic oscillator, both ^ and the presently discussed Ha,b are 
strong candidates for further investigations in this direction. 

The plan of the paper is as follows. Section 2 is devoted to describing 
elementary facts about Ha,b- In section 3 we consider a basis for which 
the eigenvalue problem associated to Ha,b is expressed as a three-term 
recurrence. The main results are to be found in section 4 where we 
establish necessary and sufficient conditions, given explicitly in terms 
of the entries of A and B, for an eigenfunction of Ha,b to be the linear 
combination of the first four elements of this basis. 



2. Elementary properties of Ha,b 

We define Ha,b rigorously as the self-adjoint operator whose domain, 
denoted below by P, is the set of all 



such that 



^,ije H\R) n |/ G L\R) : j \x^f{x)\^ < ooj = H\R)n3i]\ 
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where H'^{R) denotes the Sobolev space of index (2,2) and denotes 
Fourier transform. Since 

Ha s is a symmetric operator. It is well known (cf. [7 ) that if the 
domain of is chosen to be if^(M) fl H'^{M.), then is self-adjoint, 
non-negative and iS(M), the Schwartz space, is a core for H^- Thus 
Ha,b with domain P is a self-adjoint non-negative operator with core 
iS(R;C^). Indeed, these properties are obvious when A is a diagonal 
non-negative matrix and B = Id, the identity matrix. The general case 
follows by considering the factorization 

Ha,b = B'/'E*Hc,uEB'/^, 

where B^^^'^AB^^^'^ = E*CE is the Jordan diagonalization of the for- 
mer matrix, and by using the fact that V is invariant under the action 
of constant matrices. 

Lemma 1. The spectrum of Ha,b consists exclusively of isolated eigen- 
values of finite multiplicity whose only accumulation point is +oo. More- 
over, if \n denotes the n-th eigenvalue of this operator counting multi- 
plicity, then 

a\'%\'\2n + 1) < A2„+i < \2n+2 < al%^\2n + 1), 

where < ai < 02 and < bi < 1)2, are the eigenvalues of A and B, 
respectively. 

Proof. It reduces to showing that 

^n{HAx,Bx) ^ ^n{HA,B) < \{Ha2,B2)-, 

where Aj = aj{ld) and Bj = bj{ld). This follows directly from the 
min-max principle (cf. [8^), the estimates 

and (0). □ 
The above universal bound is not sharp in general and for most pairs 

{A,B), X2n+1 7^ A2n+2- 

As we mentioned earlier, it is not difficult to compute the eigenvalues 
and eigenfunctions of Ha,b when A and B commute. Indeed AB = BA 
if, and only if, A and B have one (and hence both) eigenvectors in 
common. Let Wj 7^ be such that Awj = ajWj and Bwj = bjWj, 
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for j = 1,2. Let be, as in (P), the eigenfunctfons of Ha- Let 

l3j = ^/ojjbj > 0. Then 

= b,{-dl + {a,/bj)x')w,(l)li^{x) 
= h]/\]/\2n + l)w,(P'',^{x). 

By choosing = 1, the family {wj(j)ri {x) : j = 1, 2; n = 0, 1, . . .} is 
an orthonormal basis of L^(]R; C^), hence 

Spec//A,B = {h]'\]'\2n + 1) : j = 1, 2; n = 0, 1, . . .}. 

The analysis below will show that finding the eigenvalues and eigen- 
functions of Ha,b whenever AB ^ BA is by no means of the trivial 
nature as the above case. 

3. Hermite expansion of the eigenfunctions in the 
non- commutative case 

Without further mention, we will often suppress the sub-indices in 
operator expressions. The structure of Ha^b = H allows us to de- 
compose L^(R;C^) into two invariant subspaces where the eigenvalue 
problem can be studied independently. We perform this decomposition 
as follows. Given a > 0, let 

ni := Span {wa;^e-°^'/2 . ^ ^ = 2k,k = 0, 1, ■ ■ ■ }, 
:= Span {t;x^e-"^'/2 : v e C\ j = 2k + 1, k = 0,1, ■ ■ ■}, 

and denote by = H\{Vnni). Since n% are invariant under d^, 
multiplication by x^ and action of constant matrices, these spaces are 
also invariant under H. Hence : V fl — > are self-adjoint 
operators and 

Spec H = Spec U Spec H~ . 

Let 

L = 2-^/\x + d^) and L* = 2-^/\x - d^) 

be the annihilation and creation operators for the scalar harmonic os- 
cillator. Then 

L^l = 0, L<p\ = n''m_, and L*cp\ = {n + if/^l^,. 
From these relations one can easily deduce the recurrent identities 

2ax^l = {n + 2fl\n + lf'^<Pl^^ + {2n + 1)^ + n"\n - if'^l^,, 
2a-'dl<Pl = {n + 2fl\n + lf'^l^,-{2n + ^W^n - 



where, here and elsewhere, any quantity with negative sub-index is 
zero. Since {0^}^o orthonormal basis for -L^(R), we can expand 
any vector of L^(R; C^) via 



(3) 



^ n=0 



for a suitable unique sequence (f„) e 1^{N; C^). Moreover, denoting by 
Na := a~^A - aB and := + aS, 



H 



[{-dl)B + x'A]J2vnrn 

n=0 



n=0 



= ^ 5^(n + 2)i/2(n + 1)1/2 jv^y^0^^^ ^ + 1)M„^;„C+ 

n=0 
k=0 

+{k + 2y/\k + iy/'N^vk+2\ ct>t- 

Thus 2H^ are, respectively, similar to the block tri-diagonal matrices 

Tt 



(4) 



S. 



± rp± 
2 -'S 



V 



acting on /^(N;C^), where 

St = {Ak + 1)M„, r+ = (2A:)V2(2A; - 1)^2 at^, 
5^- = {Ak + 3)M„ and = (2fc)i/2(2A; + l)i/2iV,. 

In order to reduce the amount of notation in our subsequent discus- 
sion, we consider Ha,b in canonical form as follows. If < 6i < 62 are 
the eigenvalues of B, let C/*(diag[6i, 62])^ be the diagonalization of B, 
and set A := h:[^UAU* and B := diagjl, 62/&i]- Then 

where 

a e 



(5) 



A 



and 



1 
b 
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Here the positivity of A and B is equivalent to the conditions 

b>l, a,c>0 and < < ac. 

Furthermore notice that AB = BA if, and only if, either 6 = 1 or 
^ = 0. Hence, unless otherwise specified, we will consider without loss 
of generality that the pair [A, B) is always the pair [A, 5) in 

By virtue of the tri-diagonal representation it seems natural to 
expect that the Hermite series Q may be a good candidate for ex- 
panding the eigenfunctions of H. Not to mention that it is the obvious 
extension of the scalar and commutative cases. In this respect, we may 
consider "finite lifetime" series expansions of eigenfunctions $ of H, 

k 

(6) <^ = J2'^nrJ'^ 

n=0 

for suitable finite /c G N U {0}, a{k) > and Vn G C^. The results we 
present below show that, contrary to the above presumption, Q is not 
such a good candidate for expanding $ for small values of k. To be 
more precise, we show that for k = 0, 1, 2, 3, an expansion of type (jU)) 
is allowed only for a small sub-manifold of the region 

(7) R := {{b, a,c,\^\) eR^ : a, c> 0, 6 > 1, < < ac} 

corresponding to all positive definite pairs {A,B). 

We first discuss the cases k = 0,1 and leave k = 2,3 for the forth- 
coming section. The following result includes a family of test bases 
larger than the one considered in ©. 

Lemma 2. Let (pni^) ^he eigenfunctions of the scalar harmonic 
oscillator Ha- Then Ha,b has an eigenf unction of the type $(a;) = 
(a(/)"(x), b4>^{x)Y where a,b E C, a, P > and m, n G N U {0}, if and 
only if AB = BA. 

Proof If if$ = A<l>, then 

-5(0°)" + ax^d^: + ixVm - = 0, 
-66(0^)" + cx^m + Ix^a^n - XHi = 0. 

If a = or 6 = in the above identities, necessarily ^ = so AB = BA. 
Hence without loss of generality we can assume that ab 0. 

If a ^ P, once again = 0. Then we may suppose that a = p. Since 
both left hand sides of the above identities are equal to p(x)e~"^ 
where in both cases p{x) is a polynomial of degree 2 -|- max(m,n). 
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necessarily either ^ = or m = n. In the latter case, the above system 
is rewritten as 

-~a{^-)"+(a~a + fb)x'^:-\~a^'^ = 0, 

-bim" + (cb + Ca)x'ra - = 0. 

Since 0^ is an eigenfunction of where a > 0, necessarily ^ G M. 
Furthermore, 

a + ^b/a = = c/b + ^a/ {bb) and 

a + ^b/a = \/{2n + l) = bc + ^ab/b. 

Hence necessarily 6=1. □ 

Since are invariant under the action of and the even (resp. 
odd) terms in the series (jHI) belong to Ti"*" (resp. 7Y~), the above lemma 
ensures that $ = fo0o + "^i^i is an eigenfunction of H if, and only if, 
A and B commute. 

4. Four-term expansion of eigenfunctions of Ha,b 

In this section we study necessary and sufficient conditions in order to 
guarantee that $ e L^(M; C^), with finite lifetime expansion of the type 
© for = 2 and 3, is an eigenfunction of H for suitable a{k) > when 
AB 7^ BA. In other words, assuming that $ satisfies the constraint 

we aim to investigate conditions ensuring = A$. 
Since 

v(j)^^ e and eHI, v e 

for all G N U {0}, and the subspaces are invariant under if, we 
may consider the even and odd cases separately. To this end, let 

(8) $+=t;o0^+t;20^ and $-=t;i0^ + t;30^, 

$^ G respectively. Then our goal is to find necessary and sufficient 
conditions, given in terms of (6, a, c, |(^|) G i? \ dR, ensuring that $^ is 
an eigenfunction of H^. The following is our main result. 

Theorem 3. Let 



-B = n L and A — , ^ 

where 6 > 1, a, c > and < < ac. Let 13 > Q he such that is 
an eigenvalue of B~^/'^AB~^/'^ > 0. Let 

_ 5(3{ab + c - 2(3%) _ 7(3{ab + c - 2(3%) 



8 L. BOULTON, S.A.M. MARCANTOGNINI, M.D. MORAN 

Then 

i) has an eigenfunction $+(a;) of type (jHI) if and only if 

(9) 2Aevcn[a + C - (6 + = 5P{P - Aeven)(Aeven - Pb) . 

In this case a = [3 and H^^^ = Aeven'^'"''- Furthermore, Aeven is an 
eigenvalue of Mjj. 

ii) H~ has an eigenfunction of type (jH)) if and only if 

(10) 6Aodd[a + c-{h+ = 7/?(3/? - Aodd)(Aodd - 3/36). 

In this case a = (3 and H^^^ = Aodd'^' ■ Furthermore Aodd is an 
eigenvalue ofiMp. 

Notice that the conditions on a, 6, c and ^ ensure that AB ^ BA. 

Proof. Put 

$+(x) = (Mo + M2x2)e-"^'/2, 

for Uo,U2 7^ 0. Since is similar to the tri-diagonal matrix (0}, then 
= A$+, if and only if 

SqUq + = Xuq 

T^Uq + S1U2 = XU2 

T+U2 = 0. 

The latter equation implies that U2 G ker and thus the first one 
implies that A is an eigenvalue of Sq = Ma with associated eigenfunc- 
tion uo- A straightforward computation shows that the above system 
is equivalent to 

(A - a^B)u2 = 

(11) (A - a^B)uo + {haB - \)u2 = 
{aB - X)uo - 2Bu2 = 0. 

The first equation holds if and only if a = (3. Here 



(12) f3 = +^ 



lab + c± v/(c-a6)2 + 4|eP6 



2b 



and U2 = { -p 1 . Notice that in this case 



i 

= det(A - p'^B) = (a - /52)(c - - 1^1'. 
'a- 

Let U2= \ 1 • Then U2 -L U2 and 



i 

{A - p^B)u2 = (a + c - (1 + b)P^)u2. 



Decompose 

UO = IU2 + 7^2, 

for suitable 7,7 G C. Then the second identity of (fTT|) holds if and 
only if, A = Aeven and 

(13) ^(a + c - (1 + h)(3^)\a - (3^) = 5/5(1 - 6)1^1^ 

The third identity of (fTT|) can be rewritten as the system 

7(/5 - A)(c - a%) + 7(/? - A) (a - (5^) = 2(c - f3%) 

7(/3&-A)(-0+7(/?&-A)(0 = -2&e 

in 7 and 7. By finding 7 from this system (for instance by Newton's 
method) and by equating to (|T!?|l . a straightforward computation yields 
Q. 

The proof of ii) is similar. □ 

Notice that there is a duality of signs in the definition of (3 (cf. fll2|) ). 
so conditions © and (fTUj) comprise two possibilities each. Let 

flf,,,,:={{b,a,c,\^\)eR : (H holds} 

and 

n^ad-={ib,a,c,\^\)ER : (HHI) holds}, 

where the sign for the super-index is chosen in concordance to the sign 
in expression ()12|) . By computing the partial derivatives of both sides of 
identities (jH)) and (jTU)), a straightforward but rather long computation 
which we omit in the present discussion, shows that these four regions 
are smooth 3-manifolds embedded in i? C M^, see ((Zj). 

The fact that f^^en ^'^^ empty is consequence of the following 
observation. By fixing 6 > 1 and |^| > 0, and by putting c = ab, 
condition Q can be rewritten as 

= b-\b + l)^[2(a6 + c - 2p%) - (/3 - A)(A - b(3)] 

= ±ffl(6+l)2-/32(l-96)(9-6), 
Vb 

where = a ± \^\b^^^'^. Then (6, a, db, |^|) G ^^^011 whenever 

. _ - m + 5) 

for 9 < 6 < 9 + 4-\/5, and (6, a, a6, |^|) G f^CTen whenever 

~ _ 1^1(562-906 + 5) 
"~ 51/2(952 _ + 9) 
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Figure 1. Left: projection of ^tven f^even onto the 
hyperplane |^| = 1, the darker color corresponds to fieven- 
Right: projection of f2^^j and ^2"^^^ onto the hyperplane 
1^1 = 1. The lighter color corresponds to 



for l<6<9or6> 9 + A^/5. Furthermore, notice that if (6, a,ab,\^\) E 
^^ven. then 

(6, ra,rc,r 1^1) G f^^en, for all r > 0. 

This shows that ^if^^^ are non-empty, unbounded and the semi-axis 
{{b, 0, 0, 0) : 6 > 1} intersects dQf^^^. All these properties also hold for 



Furthermore 



^even ^ ^odd = ^ ^ud 1],,,, H il^^^ = 0. 

Indeed, if and (fTIHl hold at the same time for the same tetrad 
(6, a, c, |,^|), then, according to TheoremEl A should be at the same time 
eigenvalue of M/3 and 3M/3 (/5 defined with the same sign). Obviously 
the latter is a contradiction. In figure 1 we reproduce the projections of 
these four regions, onto the hyper-plane |^| = 1. These picture suggest 
that nQ- ^ 0. 

Finally, we may comment on the issue of considering a more general 
basis for expanding the eigenfunctions of H. One might think that a 
natural candidate for generalizing © is the finite expansion 



4 = E 

j=0 



b,<p^ 



where Oj, bj are complex number and a, [3 > 0. We studied a particular 
case of this in Lemma El It turns out that if $ is an eigenfunction of 
the above form, then either ^ = or /3 = a, so it should be as in 
This can be easily proven by writing down the system for the eigenvalue 
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equation and considering the asymptotic behaviour of the identities as 
X — > oo. 
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